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Limits of a rapid identification of common Mediterranean sandflies 
using polymerase chain reaction-restriction fragment 
length polymorphism 
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A total of 131 phlebotomine Algerian sandflies have been processed in the present study. They belong to the spe- 
cies Phlebotomus bergeroti, Phlebotomus alexandri, Phlebotomus sergenti, Phlebotomus chabaudi, Phlebotomus 
riouxi, Phlebotomus perniciosus, Phlebotomus longicuspis, Phlebotomus perfiliewi, Phlebotomus ariasi, Phleboto- 
mus chadlii, Sergentomyia fallax, Sergentomyia minuta, Sergentomyia antennata, Sergentomyia schwetzi, Sergen- 
tomyia clydei, Sergentomyia christophersi and Grassomyia dreyfussi. They have been characterised by sequencing 
of a part of the cytochrome b (cyt b), t RNA serine and NADH1 on the one hand and of the cytochrome C oxidase 
I of the mitochondrial DMA (mtDNA) on the other hand. Our study highlights two sympatric populations within P. 
sergenti in the area of its type-locality and new haplotypes of?, perniciosus and P. longicuspis without recording 
the specimens called lex previously found in North Africa. We tried to use a polymerase chain reaction-restriction 
fragment length polymorphism method based on a combined double digestion of each marker. These method is not 
interesting to identify sandflies all over the Mediterranean Basin. 
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Algeria is a country where four leishmaniases are 
endemic. The leishmaniasis due to Leishmania infan- 
tum is transmitted by phlebotomine sandflies belonging 
to the subgenus Larroussius: Phlebotomus pernicio- 
sus, Phlebotomus perfiliewi [proven vectors according 
to Killick-Kendrick (1990)] and possibly Phlebotomus 
longicuspis (suspected vector) (Izri et al. 1990, Izri & 
Belazzoug 1993, Harrat et al. 1996, Berdjane-Brouk et 
al. 2012). Leishmania major is transmitted by the proven 
vector Phlebotomus papatasi (Izri et al. 1992). Leishma- 
nia tropica and Leishmania killicki are transmitted by 
the proven vectors Phlebotomus sergenti (Guilvard et al. 
1991, Boubidi et al. 2011, Jaouadi et al. 2012). 

The phlebotomine sandfly fauna of Algeria has been 
studied in the past (Parrot 1917, 1935, 1942, Rioux et al. 
1970a, b, Dedet et al. 1973, 1984, Dedet & Addadi 1977, 
Belazzoug & Mahzoul 1980, 1986, Belazzoug et al. 1986, 
Berchi et al. 1986, Belazzoug 1991, Russo et al. 1991). Re- 
cently, two molecular studies characterised two closely 
related species {Phlebotomus chabaudi and Phlebotomus 
riouxi) having undistinguishable or very difficultly dis- 
tinguishable females from North Africa (Bounamous et 
al. 2008, Boudabous et al. 2009) and a new species for the 
country (and for Africa) has been recorded: Phlebotomus 
mascittii Grassi (Berdjane-Brouk et al. 2011). 
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In a recent paper, Latrofa et al. (2012) suggested to 
use mitochondrial DNA (mtDNA) cytochrome b {cyt b) 
polymerase chain reaction-restriction fragment length 
polymorphism (PCR-RFLP) as a rapid molecular iden- 
tification of the common phlebotomine sandflies in the 
Mediterranean Region. The goal of this paper is to check 
the proposed method on a sampling of 17 Algerian spe- 
cies (131 specimens) belonging to three genera: Phle- 
botomus, Sergentomyia and Grassomyia. We coupled 
cyt b and cytochrome C oxidase I {COL) of the mtDNA. 
The first marker is considered as the "gold standard" 
for phlebotomine sandflies systematic. The second one 
serves as a DNA barcode for the identification of ani- 
mal species (Hebert et al. 2003). However, it was used 
in a few studies carried out on phlebotomine sandflies 
to study taxa from the Americas (Arrivillaga et al. 2002, 
Azpurua et al. 2010), from North Africa (Boudabous et 
al. 2009) and in India (Kumar et al. 2012). 

MATERIALS AND METHODS 

Sandflies collection - Sandflies were collected from 
different provinces of Algeria coupling three methods in 
order to increase the diversity (Rioux et al. 2013): sticky 
traps, CDC and ultraviolet miniature light traps and as- 
pirators (Fig. 1). They were stored in 96% ethanol. One 
hundred thirty-one males and females selected for this 
study are indicated in the Supplementary data. 

Sandflies mounting and identification - The head and 
genitalia of individual male sandflies were cut off within 
a drop of ethanol, cleared in boiling Marc -Andre solu- 
tion and mounted between slide and cover slide for spe- 
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Fig. 1: sampling realised for this study in different areas of 
Algeria. Ain-Touta is the type locality of Phlebotomus ser- 
genti. Circles: main cities; stars: locations where sandflies 
have been caught. 



cies identification. The body related to the specimen was 
stored dried in a vial at -20°C before DNA extraction. 

The specimens have been identified by observation 
of the head and genitalia under a BX50 microscope. The 
identification keys and characters used for the identifi- 
cation of specimens are those of Abonnenc (1972), Dedet 
et al. (1984), Depaquit et al. (1998a), Pesson et al. (2004) 
and Bounamous et al. (2008). Measures and photos have 
been performed using the Perfect Image software (Aries 
Company, Chatillon, France) and a video camera con- 
nected to the microscope. 

DNA extraction - Genomic DNA was extracted from 
the thorax, wings, legs and abdomen of individual sand- 
flies using the QIAmp DNA Mini Kit (Qiagen, Ger- 
many) following the manufacturer's instructions, modi- 
fied by crushing the sandfly tissues with a piston pellet 
(Treff, Switzerland) and using an elution volume of 200 
nL, as detailed in Depaquit et al. (2004). 

PCR amplification and sequencing - All the mtDNA 
amplifications were performed in a 50 \xL volume using 
5 uE of extracted DNA solution and 50 pmol of each 
of the primers. The PCR mix contained (final concen- 
trations) 10 mM Tris HC1 (pH 8.3), 1.5 mM MgCl 2 , 50 
mM KC1, 0.01% Triton X 100, 200 [iM deoxynucleotide 
triphosphate each base and 1.25 units of 5 prime Taq 
polymerase (Eppendorf, Germany). 

The cycle profiles were marker dependent. Each PCR 
begins by an initial denaturation step at 94°C for 3 min 
and finishes by a final extension at 68°C for 1 min. Am- 
plification of a fragment of cyt b gene has been done by 



using the primers N1N-PDR and C3B-PDR, following the 
method previously published by Esseghir et al. (1997): 
five cycles (denaturation at 94°C for 30 s, annealing at 
40°C for 60 s and extension at 68°C for 60 s) followed by 
35 cycles (denaturation at 94°C for 60 s, annealing at 44°C 
for 60 s and extension at 68°C for 60 s). Their COI domain 
was amplified using the primers used by Hajibabaei et al. 
(2006): LepF and LepR, under the following thermal pro- 
file (Costa et al. 2007): five cycles (denaturation at 94°C 
for 30 s, annealing at 45°C for 90 s and extension at 68°C 
for 60 s), then 35 cycles (denaturation at 94°C for 30 s, 
annealing at 51°C for 90 s and extension at 68°C for 60 s). 

Amplicons were analysed by electrophoresis in 1.5% 
agarose gel containing ethidium bromide. Direct sequenc- 
ing in both directions was performed by Sanger's method 
using the primers used for DNA amplification. Reagents 
for PCR cleanup were Agencourt® AMPure XP-PCR 
Purification and those for sequencing reaction were Big 
Dye Sequencing Buffer (Applied Biosystems Foster City, 
California, USA); Agencourt® CleanSEQ kit: CleanSeq® 
beads. The instruments used for the sequencing were 
Biomek NXp BiomekFXP (Bekman Coulter); Sequencer 
3730 XLT (Applied Biosystems). The correction of se- 
quences was done using Pregap and Gap softwares in- 
cluded in the Staden Package (Bonfield & Staden 1996). 

Molecular analyses - They are based on the two 
datasets of sequences. Analyses were performed using 
haplotypes obtained from this study and sequences of P. 
chabaudi and P. riouxi available in GenBank (Supple- 
mentary data). Sequence alignment was performed us- 
ing the CLUSTALW routine included in the MEGA v.4 
software (Tamura et al. 2007) and checked by eye. The 
objective of this study is to provide a database for sand- 
flies identification, not to do any phylogenetical analy- 
sis. Consequently, a neighbour-joining (NJ) analysis was 
performed using MEGA 4 software, with the Kimura-2 
parameters model and using uniform rates among sites. 
Gaps were treated as missing data. 

RFLP - The diagnostic endonuclease restriction sites 
on cyt b and COI mtDNA sequences were predicted 
for each specimen using CLC workbench 5.2 software 
(clc-genomics-workbench.com-about.com). A panel of 
restriction enzymes was tested including Asel enzyme 
proposed by Latrofa et al. (2012) for cyt b digestion. Be- 
cause one restriction enzyme cannot provide an original 
digestion pattern per species, we selected double diges- 
tion for both molecular markers. 

PCR-RFLP assays were performed in a 50-uE total 
volume reaction mix, containing 15 uL of PCR product 
(from PCR vials), 0.1 uL of Asel, 0.4 |iL of Mnll for cytb 
and 0.05 uE of Mspl, 0.05 p.L of Taql for COI, 5 uL of 
NEB buffer 3 for cyt b and 5 uL of NEB buffer 4 for 
COI containing bovine serum albumin (New England 
Biolabs, France). 

For cyt b PCR-RFLP, we selected a double digestion 
coupling Asel with Mnll. PCR products were digested 
during 2 h at 37°C. 

For CO/PCR-RFLP, we selected a double digestion cou- 
pling Mspl with Taql. According to their different tempera- 
ture of activity (37°C and 65°C, respectively), PCR products 
were digested during 1 h at 37°C then 1 h at 65°C. 
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The digested samples were separated by electropho- 
resis in a 3% agarose gel to produce DNA fragments and 
sized by comparison with markers 50 bp ladder, 100 bp 
ladder and 20 bp ladder (Cliniscience, France). 

According to the sequencing of all the specimens 
processed in the present study, we did not use any posi- 
tive control, but we checked that each restriction enzyme 
had functioned properly in each reaction by comparison 
to the predicted digestion. 

RESULTS 

PCR amplification was successful for all the speci- 
mens processed. The GenBank accessions for COI 
and cyt b are indicated in the Supplementary data. The 
length of the analysed markers is of 680 bp for COI. It 
varies from 510-525 bp for cyt b. Each gene from each 
haplotype had an open reading frame (ORF). The se- 
quences labelled COI include exclusively this marker. 
The sequences labelled "cytB" included in fact an ORF 
of cytB (positions 1-321), the t RNA serine (positions 
321-378), then the ORF for NADH subunit 1 (from 379 to 
the last position). These ORFs are translated in proteins, 
explaining the low probability they could be pseudo- 
genes (Rogers & Griffiths-Jones 2012). 

Global trees based on cyt b and COI sequences are 
presented in Figs 2, 3, respectively. According to De- 
paquit et al. (1998b), Phlebotomus bergeroti has been 
selected to root the tree. All the species morphological- 
ly recognised are well individualised. We note that P. 
sergenti included two populations without any morpho- 
logical difference. The identification of P. perniciosus 
and P. longicuspis is not doubtful. We did not record 
any atypical specimen. 

The double digestion of each PCR product confirms 
the expected fragments for all haplotypes (Tables I, II). 
The resolution of DNA fragment size by the gel fraction- 
ation method used is about 10 bp. The double digestion 
of COI by restriction enzymes TaqI and Mspl is very ef- 
ficient, but cannot unequivocally distinguish Phleboto- 
mus ariasi and Sergentomyia schwetzi. 

The double digestion of cyt b by restriction enzymes 
Asel andMnll provide restriction profiles indicated in Table 
II. On the one hand, it cannot individualise one population 
of P. sergenti, P. ariasi and P. riouxi. On the other hand, it 
cannot distinguish P. longicuspis and P. perfiliewi. 

No partial digests were recognised in the analysis. 

DISCUSSION 

The specimens identified by morphology as belong- 
ing to a species are branched together regarding inde- 
pendently cyt b or COI sequences (Figs 2, 3). 

Cyt b sequences provide a NJ tree in agreement with 
the traditional morphological taxonomy of the phle- 
botomine sandflies (Fig. 2). Sergentomyia are grouped 
together including Grassomyia dreyfussi. Concerning 
Larroussius, all of them are grouped together and two 
branches are individualised: one including P. ariasi and 
Phlebotomus chadlii and another containing P. perni- 
ciosus, P. longicuspis (including 2 lineages) and P. per- 
filiewi. These data are in accordance with those obtained 
by Esseghir et al. (2000) on cyt b and Di Muccio et al. 



(2000) on rDNA internal transcribed spacer 2. Concern- 
ing the subgenus Paraphlebotomus, the species P. ser- 
genti, P. chabaudi and P. riouxi are grouped together. 
The first one shows two lineages in its type-locality. 
Moreover, Phlebotomus alexandri is not included in this 
branch, as previously observed (Depaquit et al. 2000, 
Kriiger etal. 2011). 

The NJ tree based on COI sequences (Fig. 3) has a 
surprising topology: the species belonging to the sub- 
genera Sergentomyia, Larroussius and Paraphlebotomus 
are not grouped together. Despite this curious branching, 
some results are congruent with cyt b: (i) the position of 
P. alexandri, (ii) the existence of two molecular lineages 
within P. sergenti topotypes, (iii) the individualisation of 
P. chabaudi and P. riouxi, (iv) the existence of two lineag- 
es within P. longicuspis and (v) a high variability within 
Sergentomyia minuta and Sergentomyia antennata. 
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Fig. 2: neighbour-joining tree based on mitochondrial DNA cyto- 
chrome b sequences rooted on Phlebotomus bergeroti. Bootstrap 
values after 1,000 replicates are indicated on the branches. 
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Fig. 3: neighbour-joining tree based on mitochondrial DNA 
cytochrome C oxidase I sequences rooted on Phlebotomus 
bergeroti. Bootstrap values after 1,000 replicates are indicated 
on the branches. 
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Many lineages have been identified in P. sergenti in 
the literature from populations from different parts of 
the species distribution area and no study emphasise a 
link between the molecular variability and the morphol- 
ogy (Depaquit et al. 2002, Yahia et al. 2004, Moin-Vaziri 
et al. 2007, Baron et al. 2008, Dvorak et al. 2011). The 
two mitochondrial lineages {cyt b as well as COI) within 
Algerian specimens of P. sergenti coming from different 
localities, all located just around the type locality area 
(Fig. 1) called Ain Touta, formerly Mac Mahon (Parrot 
1917). This locality has not been precisely designated by 
Parrot (1917). The specimens processed in the present 
study can be considered as being topotypes, clearly la- 
belled and stored in the collection of the laboratory of 
Parasitology of the Faculty of Pharmacy of Reims. These 
two populations are strongly separated and are character- 
ised by many variable nucleotidic positions: about 30 for 
cyt b and 40 for COI (Fig. 4). The mean pairwise distance 
between the two populations of P. sergenti (> 5%) is com- 
parable to the pairwise distances individualising P. perni- 
ciosus from P. longicuspis or P. perfiliewi in the present 
study and question about the status of these populations. 

The haplotypes obtained for P. perniciosus and P. 
longicuspis have been compared with those of cyt b 
available in the literature (Esseghir et al. 1997, Pesson 
et al. 2004, Perrotey et al. 2005) (Fig. 5). The main P. 
perniciosus Algerian haplotype is the same than the 
main Mediterranean haplotype (= pernOl), but three new 
haplotypes are recorded from Algeria for this species. 



Concerning P. longicuspis, the most common haplotype 
is the lcusOl. However, we record four new haplotypes 
including a couple (LC2 and LC637) strongly individu- 
alised from the other ones. The COI sequences also in- 
dividualise the latter specimens from all other P. longi- 
cuspis (Fig. 3). 

Recently, the single digestion based on cyt b sequenc- 
es by Latrofa et al. (2012) as a rapid molecular identifi- 
cation method for the common phlebotomine sandflies 
in the Mediterranean Region does not apply in Algeria. 
In fact, these authors focused on five species only, com- 
monly caught in Italy: P. papatasi, P. perniciosus, P. per- 
filiewi, Phlebotomus neglectus and S. minuta. We tried 
this method on the Algerian sandflies, including P. pa- 
patasi. The method is not able to distinguish all the spe- 
cies and the combined double digestion of two different 
markers is needed for the specific identification. Due to 
the conservation of some parts of DNA sequences, the 
simple digestion of PCR products cannot separate some 
species belonging to different genera (like P. ariasi and 
S. schwetzi) or subgenera (like P. sergenti and P. ariasi). 
Consequently, this method is not enough discriminant to 
be used in routine all over the Mediterranean Basin. In 
fact, it is easier to firstly identify the species easy to rec- 
ognise by a microscopical examination and secondly, to 
apply efficient PCR-RFLP methods to identify the species 
for which morphological identification is difficult, like the 
females of the Perniciosus complex, those of P. chabaudi 
and P. riouxi and for some Sergentomyia males. 



TABLE I 

Species and populations examined in the present study - number and positions of cuts predicted 
for cytochrome C oxidase I mitochondrial DNA restriction fragment length polymorphism 



Taql and Msp\ 
(expected) 





Cut 


Fragments 




Species 


(n) 


(n) 


Cut sites 


Phlebotomus bergeroti 


4 


5 


85,91, 101, 176, 254 


Phlebotomus alexandri 


5 


6 


22, 84, 85, 91, 156, 267 


Phlebotomus longicuspis 


4 


5 


22,73, 91, 113,408 


Phlebotomus perniciosus 


3 


4 


22,91, 184, 408 


Phlebotomus chabaudi 


5 


6 


30, 85, 89, 91, 158, 256 


Phlebotomus sergenti (group 1) 


3 


4 


25, 62, 111, 507 


Phlebotomus sergenti (group 2) 


3 


4 


22, 85, 91, 507 


Phlebotomus chadlii 


6 


7 


14, 22, 85, 91, 101, 174, 224 


Phlebotomus riouxi 


6 


7 


21,30, 85,91, 141, 157, 186 


Phlebotomus ariasi 


2 


3 


85,91,527 


Phlebotomus perfiliewi 


2 


3 


11, 184, 408 


Sergentomyia fallax 


2 


3 


85, 265, 353 


Sergentomyia schwetzi 


2 


3 


85,91,527 


Sergentomyia clydei 


4 


5 


82,91, 174, 176,184 


Sergentomyia minuta 


3 


4 


88, 90, 176, 350 


Sergentomyia christophersi 


5 


6 


32,81,85,91, 101,319 


Sergentomyia dreyfussi 


3 


4 


70, 91, 184, 360 


Sergentomyia antennata 




4 


82, 110, 157, 356 
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TABLE II 

Species and populations examined in the present study - number and positions of cuts predicted 
for cytochrome b mitochondrial DNA restriction fragment length polymorphism 



Asel and Mnll (expected) 





Cut 


Fragments 




Species 


(n) 


(n) 


Cut sites 


Phlebotomus bergeroti 


3 


4 


14, 20, 74, 434 


Phlebotomus alexandri 


5 


6 


1 1 "\ J \ S -\ |" \ "\ 1/1(1 "V ,1 j 1 

14, 20, 62, 92, 108, 248 


Phlebotomus longicuspis 


5 


6 


14, 20, 62, 92, 103, 263 


Phlebotomus perniciosus 


4 


5 


14, 20, 62, 103, 353 


Phlebotomus chabaudi 


4 


5 


14,20, 62, 198,248 


ni ill * • f 1 \ 

Phlebotomus sergenti (group l) 


4 


5 


14, 20, 96, 208,209 


Phlebotomus sergenti (group 2) 


3 


4 


14, 20, 62, 449 


Phlebotomus chadlii 


3 


4 


20, 62, 107, 356 


Phlebotomus riouxi 


3 


4 


14,20, 62, 445 


Phlebotomus ariasi 


3 


4 


14, 20, 62, 453 


Phlebotomus perfiliewi 


5 


6 


14, 20, 62, 92, 103, 263 


Sergentomyia fallax 


3 


4 


20, 62, 103, 353 


Sergentomyia schwetzi 


4 


5 


20, 52, 52, 62, 356 


Sergentomyia clydei 


5 


6 


10, 20, 62, 103, 134,212 


Sergentomyia mimita 


6 


7 


14, 19,20, 42,74, 93,296 


Sergentomyia christophersi 


6 


7 


20, 34, 46, 58, 62, 134, 184 


Sergentomyia dreyfussi 


7 


8 


14, 15, 20, 28, 62, 103, 134,171 


Sergentomyia antennata 


1 


2 


20, 515 



Cytochrome B: 

< CytE >< tRHA-Ser >< NAD HI > 
12222222223333333333333444444 

457999313157768901133 5 67789992 55566 

27878981279470 6773 6142329812 4714703 
SE304 AGGATATATAGTGGACGTGGGTTAGGTTAAGGTCC 

SE307 

SE312 C....A....C 

SZ335 C 

5E669 C....A....C A.. 

SE670 C....A....C 

SE675 A....CC 

SE305 GAAGGGCGAGA.AAGTACAAA. .GAAC.GTAA.TT 
SE394 GAAGGGCGAGA.AAGTACAAA. . GAACCGTAA. TT 
5E539 GAAGGGCGAGA . AAGTACAAA . . GAACCGTAA. TT 
SES40 GAAGGGCGAGA.AAGTACAAA. .GAAC.GTAA.TT 
SE621 GAAGGGCGAGA.AAGTACAAA. .GAAC.GTAA.TT 
SE627 GAAGGGCGAGA.AAGTACAAA. .GAAC.GTAA.TT 

COL 

< coi > 

11111222222222333333333334444445SSSSS5S56666 
5678899028890 133345790 133 44778991 12 3 4801334SSS8122S 
35989473547 627258708840 173969S14S873 840S23714701282 
SE304-COI ACCATTTAATTTTTGTCACCTATAGAATTTACCTATAACTAACTTCGTTCG 

5E307-COI G 

SE312-COI G 

SE335-COI . A 

SE669-COI . . . G G G 

SE67Q-COI G G A 

SE675-COI G G G. A 

SE3Q5-COI GTT. CCC. . . C . CCACTCTTC. C . A. GCCCCTTCGCG. TCGGTC . TACCTT 
SE394-COI GTT . CCCGG . C . CCACTCTTC . G . AGGCCCCTTCGCG. TCGGTC . TACCTT 
SE539-COI GTT . CCCGG . C . CCACTCTTC . G. AGGCCCCTTCGCG. TCGGTC . TACCTT 
SE540-COI GTT.CCCG. . C . CCACTCTTC . C . A. GCCCCTTCGCG. TCGGTC . TACCTT 
SE627-COI GTT. CC . G . . C . CCACTCTTC . C . A. GCCCCTTC . CG. TCGGTC . TACCTT 
SE621-COI GTT . CCCG . ACACCACTCTTC . C . A. GCCCCTTCGCGGTCGGTCCTACCTT 
*** ** * ***** *** * ******** ** *** + * ***** 

Fig. 4: variable positions observed within Phlebotomus ser- 
genti for cytochrome b and cytochrome C oxidase I (COI) mi- 
tochondrial DNA. Stars indicate the sites characterising the 
two populations. 



pernOl (France, Italy, Malta, Tunisia, Morocco} 
pern[>2 (Italy) 
pernCH (Spain) 
pern05 (Spain) 

parnOMPNA: Tunisia, Morocco) 
pern07(PWA: Morocco) 
pernOS (France) 

PN-main-algerian-haplotype (n-16) 

PNFE1 . 2-CYT 

PN592-CYT 

PNG42-PN643 

lcusOl (Tunisia, Morocco) 
Lous 02 (Morocco) 

LC3-549-601-BE1-54 5-614-509-54 6-553-331 

LC633 

LC2 

LC637 

LCVIM6 

IcxOl (Morocco) 



11111111122222 
124 55777 734 4 6678 &90 14 55 
67 9579 4 824 58357 58113714 628 
TAAAATTTTCTCTftC/ATCGACfiGATG 



.G. TCCCC. .T.GT.CT.C. . . TCA 
. G . TCCCC . . T . GT . CTAC . . . TCA 
.G. TCCCC. .T.GT.CT.C. . . TCA 
. G . TCCCC . .T|GT . CT . C . . . TCft 

.g.tcIcc. .t.g|.c1.c.|.tca 

.G.TclcC. + T.G|.CHC.|.TCft 
. G . Tc|cc . . T . GT . CT . C . . . TCA 

. CC.CTC. . .Oh T.fi 



Fig. 5: alignment of cytochrome b haplotypes concerning 
Phlebotomus perniciosus and Phlebotomus longicuspis and 
showing exclusively the variable sites reported by Esseghir et 
al. (1997), Pesson et al (2004) and Perrotey et al. (2005). The 
Algerian haplotypes are written in bold. The variations ob- 
served in Algerian specimens are underlined. The stars indi- 
cate the sites individualising the haplotypes LC2 and LC637. 
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